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ABSTRACT 

A  test  program  to  evaluate  a  new  transonic  axial 
compressor  stage  was  conducted.  The  stage  was  designed  (by 
Nelson  Sanger  of  NASA  Lewis)  relying  heavily  on  CFD 
techniques  while  minimizing  conventional  empirical  design 
methods.  The  stage  was  installed  in  the  NFS  Transonic 
Compressor  Test  Rig  and  instrumented  with  fixed  temperature 
and  pressure  probes.  A  new  PC-based  data  acquisition  system 
was  commissioned  and  programmed  for  stage  performance 
measurements.  These  were  obtained  at  50,  60,  65,  70,  and  80% 
of  the  design  speed  before  failure  of  the  spinner  retaining 
bolt  led  to  the  loss  of  the  stage.  The  flow  through  the 
rotor  was  analyzed  and  the  rotor  perfozmnance  predicted  using 
a  3-dimensional  viscous  code  (RVC3D) .  The  predicted  rotor 
performance  agreed  qualitatively  and  was  niimerically 
consistent  with  the  measured  stage  performance. 
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I .  INTRODUCTION 


A  new  transonic  axial  compressor  stage  designed 
specifically  for  the  Naval  Postgraduate  School  (NPS) 
Transonic  Coitpressor  Rig  (TCR)  was  conpleted  at  NASA  Lewis 
in  1994,  by  Nelson  L.  Sanger  [Ref.  1].  A  schematic  of  the 
TCR  is  shown  in  Fig.  1.  The  new  stage  design  was  developed 
to  serve  as  a  vehicle  for  instructional  purposes,  and  to 
provide  a  test  bed  for  research.  With  the  emergence  of 
Conputational  Fluid  Dynamics  (CFD)  as  a  reliable  tool,  the 
design  of  the  stage  developed  into,  as  Sanger  described  it, 
"...the  ultimate  CFD  validation  experiment".  The  design 
parameters  produced  from  this  process  are  presented  in  Table 
1. 

The  manufactured  stage,  after  assembly  in  the  TCR,  is 
shown  in  Fig.  2.  The  stage  was  installed  in  the  same  manner 
as  previously  operated  stages  and  as  documented  in  Ref .  2 . 
However,  a  Lucite  case  wall  was  installed,  as  seen  in  Fig. 
3,  and  a  new  data  acquisition  system  was  enployed.  The  data 
acquisition  system  used  a  Personal  Conputer  (PC)  running 
Hewlett-Packard  HPVEE  software  to  control  various 
instruments  in  a  HP75000  VXI-bus  mainframe.  The  stage  was 
tested  and  measurements  taken  at'  50,  60,  65,  70  cind  80 

percent  of  design  speed.  Testing  was  unavoidably  halted 
during  the  test  at  80%  speed  due  to  the  failure  of  the 
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spinner  retaining  bolt,  which  resulted  in  the  loss  of  the 
stage. 

The  present  report  documents  the  first  test  program  of 
the  Sanger  stage.  The  report  contains  three  main  sections. 
The  first  section  details  the  rig  and  the  newly  developed 
data  acquisition  system.  The  second  section  focuses  on  the 
results.  The  final  section  gives  the  results  of  a  CFD 
analysis  of  the  rotor  using  the  "Rotor  Viscous  Code  3 
Dimensional"  (RVC3D)  code  developed  by  R.  Chima  [Ref.  3]. 


Parameter 

Final 

Selection 

Rotor  Pressure  Ratio 

1.61 

Stage  Pressure  Ratio 

1.56 

Tip  Speed 

1300  ft/sec 

Design  Mass  Flow 

17.09  Ibm/sec 

Specific  Mass  Flow 

Specific  Head  Rise 

0.246 

Tip  Inlet  Relative  Mach 

Number 

1.28 

Hub /Tip  Radius  Ratio 

0.51 

Rotor  Inlet  Ramp  Angle 

28.2  degrees 

457  HP 

Table  1.  Design  Parameters  for  the  Sanger  Stage. 
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Fig.  1  Schematic  of  the  Transonic  Coirpressor  Test  Rig. 
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II.  TRANSONIC  COMPRESSOR  RIG  AND  INSTRUMENTATION 


A.  COMPRESSOR  INSTALLATION 

The  general  layout  of  the  TCR  is  shown  in  Fig.  1.  The 
manner  in  which  the  stage  was  mounted  on  the  rig  was 
unchanged  from  the  previously  installed  Vavra  stage 
documented  in  Ref.  2.  The  flow  entered  the  rotor  axially.  A 
honeycomb  section  ensured  that  the  flow  left  the  cortpressor 
with  zero  circiimferential  velocity.  In  this  manner  the 
torque  supplied  to  the  rotor  was  equal  to  the  torque 
experienced  by  the  blading  and  honeycomb  of  the  stator 
section.  The  stator  section  was  mounted  on  ball  races  and 
was  free  to  rotate  against  a  flexure  instrumented  with 
strain  gauges.  The  strain  gauge  reading  was  then  a  measure 
Of  the  torque  supplied  to  the  rotor. 

The  tip  clearance  between  the  rotor  and  the  steel 
casing  that  had  been  used  previously  was  a  concern.  The 
decision  was  made  to  use  a  similar  clear  plastic  (Lucite) 
casing  until  the  steel  casing  could  be  modified  to  have  an 
abradable  material  where  the  rotor  might  anab.  The  Lucite 
casing  proved  to  fxmction  extremely  well.  The  static  tip  gap 
clearance  was  0.010-0.012  inches  which  was  larger  than  the 
0.004  inches  the  design  called  for;  however,  this  ensured 
that  no  rxob  would  occur  during  the  first  test  of  this  unique 
stage.  Additionally,  the  transparent  case  provided  an 
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opportiinity  to  see  transonic  turbomachinery  in  operation. 
One  interesting  observation  was  that  excess  oil  from  the 
bearing  lubrication  system  indicated  a  flow  pattern  and  was 
seen  to  collect  between  the  casing  and  the  rotor  tip. 

B.  INSTRUMENTATION 

1 .  Measurement  Devices 

Fig.  4  shows  the  location  of  measurement  stations 
relative  to  the  stage.  The  instrimientation  for  stage 
performance  measurements  included  Kiel  pressure  probes, 
Kiel /thermocouple  combination  probes,  static  pressure  taps, 
torque  flexure,  flow  nozzle,  and  a  magnetic  pick  up.  Total 
pressure  was  measured  using  United  Sensor  Kiel  Probes.  Two 
types  were  used,  a  1/16  inch  "Miniature  Head"  Kiel  Probe 
(United  Sensor  KAA-8,  Ref.  4)  and  1/8  inch  "Standard  Head" 
combination  Kiel /thermocouple  Probe  (United  Sensor  KT-8-J- 
12-C,  Ref.  4)  .  The  1/16  inch  mini  probe  had  a  yaw  range  of 
52°  and  pitch  range  of  40°.  Range  was  defined  by  the 
manufacturer  as  the  point  at  which  error  was  equal  to  1%  of 
dynamic  pressure.  The  1/8  inch  standard  probe  had  a  yaw 
range  of  48°  and  a  pitch  range  45°.  The  static  pressure  ports 
were  1/32  inch  taps  located  in  the  wall  along  the  length  of 
the  casing,  and  in  the  hvib  of  the  stator  at  stations  2  and 
3. 
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Thermocouples  were  J  type  ( iron-const  ant  an)  in  a  probe 
combined  with  a  Kiel  type  in5)act  probe.  The  thermocouple  had 
an  exposed  junction  and  bleed  slot. 
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Fis  •  4  Schematic  of  the  Flow  Path  Showing  Measurement 

Stations . 


The  torque  acting  on  the  stator  and  exit -honeycomb 
section  (equal  to  the  rotor  torque)  was  measured  by  the 
deflection  of  a  cantilevered  flexure.  The  cantilever  arm  was 
instrumented  with  four  strain  gauges  providing  a  near  linear 
response  to  bending.  The  arrangement  and  calibration  method 
of  the  torque  measurement  system  is  extensively  documented 
in  Ref.  2.  Calibration  of  the  torque  measuring  system  was 
necessary  each  time  the  rig  was  rebuilt.  The  calibration 
cuirve  obtained  in  the  present  study  and  resulting  equation 
used  in  the  data  reduction  program  are  shown  in  Fig.  5. 

Rotational  speed  was  measured  using  a  magnetic  flux 
cutter  and  a  disk  with  six  holes  moiinted  on  the  shaft  of  the 
drive  turbine  (Fig.  1) .  This  device  produced  six  pulses  per 
revolution  making  the  frequency  (in  Hz)  from  the  pickup 
multiplied  by  10  equal  to  the  revolutions  per  minute (RPM)  of 
the  stage. 

Mass  flow  rate  was  measured  using  a  flow  nozzle  located 
upstream  of  the  stage  as  shown  in  Fig.  1.  Mass  flow  rate  was 
calculated  using  the  pressure  difference  of  the  total 
pressure  (p^)  and  static  pressure  (p^)  at  the  flow  nozzle  and 
the  total  temperature  of  the  flow.  The  differential  pressure 
between  p^  and  p^  was  displayed  on  a  water  manometer  in  the 
control  room  to  allow  precise  control  of  mass  flow  rate 
settings  during  the  throttling  process. 
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2 .  Instrument  Placement 

Two  combination  probes  were  placed  upstream  of  the 
rotor,  each  two  inches  from  the  wall  of  the  casing  to 
measure  inlet  total  pressure  and  temperature.  At  station  3, 
17  Kiel  probes  and  three  Kiel/  thermocouple  combination 


probes  were  placed  as  shown  in  Table  2 . 


Probe 

# 

Radial  Distance 

(y) 

Angle 

18 

Kiel 

0.7 

15 

15 

Kiel 

0.4 

30 

12 

Kiel 

0.1 

45 

41 

Static  Port 

n/a 

60 

11 

Kiel 

0.05 

75 

29 

Kiel 

1.8 

90 

■ 

Positive  \ 

Angle  \ 

26 

Kiel 

1.5 

105 

23 

Combo 

1.2 

120 

20 

Kiel 

0.9 

135 

17 

Kiel 

0.6 

150 

14 

Kiel 

0.3 

165 

42 

Static  Port 

n/a 

180 

/ 

■n 

n/a 

n/a 

195 

\  / 

28 

Kiel 

1.7 

210 

25 

Kiel 

1.4 

225 

22 

Combo 

1.1 

240 

19 

Kiel 

0.8 

255 

16 

Kiel 

0.5 

270 

13 

Kiel 

0.2 

285 

43 

Static  Port 

n/a 

300 

30 

Kiel 

1.85 

315 

27 

Kiel 

1.6 

330 

24 

Kiel 

1.3 

345 

21 

Combo 

1 

360 

Table  2 ,  Location  of  Probes  at  the  Stage  Exit . 
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The  pressures  acquired  from  the  inpact  probes  at 
station  3  (probe  numbers  11-30  in  Table  2)  were  mass 
averaged  to  obtain  the  total  pressure  at  the  exit  of  the 
stage. 

Static  pressure  ports  were  placed  at  stations  Inf,  1,  2 
and  3  along  the  Lucite  casing  (Fig.  4) .  Additionally,  static 
taps  were  drilled  in  the  stator  to  measure  hiib  static 
pressure  at  station  2.  The  exit  spacer  between  the  stator 
and  honeycomb  also  had  four  static  taps  to  measure  exit 
static  pressure  on  the  hiib  at  station  3 . 

C.  TEST  PROCEDURE 

The  test  compressor  was  driven  by  two  opposed-rotor  air 
turbine  stages  (on  a  common  shaft)  supplied  from  the 
laboratory  conpressed-air  power  supply.  The  laboratory 
coirpressed-air  power  supply  was  a  12-stage  All  is -Chalmers 
axial  conpressor  supplying  up  to  11  pounds  per  second  at  up 
to  two  atmospheres  gauge  pressure.  The  test  conpressor  was 
connected  by  a  splined  shaft  to  the  rotor  of  the  turbine 
drive  unit.  The  speed  of  the  drive  \init  was  controlled  using 
an  electrically  actuated  butterfly  valve.  Air  was  drawn  into 
the  test  conpressor  through  a  rotating  plate  throttle 
followed  by  a  settling  chamber  containing  perforated  plates 
(Fig.  1).  A  flow  nozzle  was  positioned  upstream  of  the 
rotor.  The  calibration  of  the  flow  nozzle  was  documented  in 
Ref .  5 . 
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Preparation  for  a  run  consisted  of  bringing  the  Allis- 
Chalmers  coitpressor  up  to  stable  speed,  verifying  that 
conpressed  air  was  available  to  the  balance  piston  (which 
maintained  low  axial  loading  on  the  coiipressor  bearings) , 
and  setting  the  lubrication  system.  The  bearing  lubrication 
system  was  an  oil-mist  system.  The  operating  tenperature  of 
the  bearings  was  found  to  be  directly  related  to  the  supply 
pressure  of  the  air  driving  the  oil  mist.  Increasing  the 
pressure  of  the  oil  mist  supply  air  had  the  observed  effect 
of  lowering  the  bearing  operating  tenperature.  Adjusting  the 
oil  droplet  rate  down  to  2-3  drops -per-minute  was  essential. 
According  to  the  bearing  manufacturer,  wetting  the  bearing 
well  before  a  run  was  of  extreme  importance.  Allowing  the 
oil  mist  to  operate  for  fifteen  minutes  before  operating  was 
the  procedure  adopted. 

During  a  compressor  test,  operating  data  were  monitored 
by  the  operator  at  the  control  console.  Records  were  kept  of 
oil  drop  rate  and  pressure  to  the  bearing  ,  lubrication 
system,  of  the  tenperatures  of  the  turbine  and  conpressor 
bearings,  and  of  vibration  levels  from  four  accelerometers 
on  the  turbine  drive  \init.  The  axial  force  was  constantly 
monitored  and  the  balance -pis  ton  air  adjusted  to  keep  the 
axial  force  constant  at  a  level  approximately  equal  to  the 
force  produced  at  5000  RPM.  These  data  were  mainly  of 
interest  for  the  safety  and  operation  of  the  unit  and  will 
not  be  reported. 
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A  typical  test  would  be  conducted  at  a  constant  speed. 
The  run  started  with  the  coitpressor  throttle  wide  open. 
After  the  speed  had  settled  to  the  desired  level,  a 
measurement  recording  cycle  was  started.  Then  the  throttle 
would  be  closed  fractionally  imtil  the  flow  was  reduced  to 
the  desired  amount.  Then  the  drive  turbine  speed  was 
adjusted  until  it  was  steady  again  and  at  the  desired  level. 
This  procedure  was  repeated  iintil  stall  occurred.  After  the 
stall  point  was  reached  (up  to  70%  speed  only)  the  process 
was  progressively  reversed  to  obtain  intermediate  settings 
until  the  throttle  was  fully  open. 
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III.  DATA  ACQUISITION  AND  REDUCTION 


A.  OVERVIEW 

In  order  to  standardize  equipment  for  data  acquisition 
in  the  three  buildings  at  the  NPS  Turbopropulsion  Laboratoiry 
the  decision  was  made  to  use  the  HP75000  Series  B  VXI-Bus 
Mainframe  controlled  by  HPVEE  software  running  on  a  personal 
coitputer  (PC)  .  This  allowed  the  Hewlett-Packard  (HP) 
instruments  with  the  HP-IB  interface  that  had  been  used 
previously  to  remain  viable  while  incorporating  a  new  state- 
of-the-art  data  acquisition  system.  In  the  present  work  the 
mainframe  with  an  internal  digital  voltmeter  (DVM) ,  two 
scanning  multiplexers,  a  switchbox  multiplexer,  a  co\inter, 
and  a  Quad  8-bit  Digital  I/O  Module  were  used  exclusively. 
The  system,  shown  in  Fig.  6,  proved  to  be  flexible,  easy  to 
interrogate  when  connecting  new  sensors,  and  provided  real 
time  feedback  when  acquiring  data. 

B.  HARDWARE 

The  use  of  the  HP75000  Mainframe  allowed  a  variety  of 
instruments  to  be  grouped  together  and  addressed  directly. 
Communication  between  the  PC  (with  a  HP  82341C  Controller 
Card  installed)  and  the  mainfrcime  was  via  a  HP-IB  (IEEE-488) 
interface  cable.  The  fxinctions  performed  with  the 
instruments  in  the  mainframe  were  Scanivalve  control, 
measurement  of  Scanivalve  transducer  output,  thermocouple 
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measurement , 
coianting. 


strain 


gauge  measurement 


and  frequency 


Fig.  6  Schematic  of  the  Data  Acquisition  System. 
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1 .  Scan! valve  Control 

The  Scan! valve  and  Scanivalve  control  (HG-78)  system 
that  the  present  system  was  designed  to  replace,  are 
described  in  Ref.  6.  Scanivalve  control  involved  stepping 
and  homing  the  48-port  pneumatic  valve  and  ensuring  that  the 
correct  port  was  selected  and  measured.  One  HP  E1345A  16- 
channel  multiplexer  was  configured  as  a  switchbox  as  shown 
in  Fig.  7.  This  configuration  allowed  the  control  of  four 
different  Scanivalves  with  one  switchbox.  The  measurements 
reqpiired  for  the  test  program  required  the  control  of  only 
one  Scanivalve  at  position  1.  To  reset  the  Scanivalve, 
channel  12  was  closed,  sending  the  signal  to  the  HG-78  to 
reset  the  Scanivalve  connected  to  position  1.  To  step  to  the 
next  port,  channel  8  was  closed,  sending  the  signal  to 
advance  one  port.  The  position  of  the  Scanivalve  was 
obtained  from  the  HG-78,  which  contained  LED  displays  of  the 
port  setting,  manual  step  and  home  controls  (in  parallel), 
and  signal  conditioning  for  the  Scanivalve  transducers.  From 
the  HG-78,  eight  TTL  lines  (from  each  Scanivalve)  were 
connected  to  one  of  four  addressable  sections  of  the  HP 
E1330B  Quad  8-bit  Digital  I/O  card.  The  port  which  the 
Scanivalve  was  currently  set  at  was  then  read  (using 
negative  logic)  and  converted  to  a  decimal  nximber  using 
Binary  Convert  Decimal  (BCD)  logic.  The  assignment  of  the  48 
Scanivalve  ports  is  shown  in  Table  3 . 
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Port# 

Scanivalve  1 

1 

Tare 

2 

Scale 

3 

Flow  Nozzle  P6 

4 

Flow  Nozzle  Ps 

5 

Inlet  Ptinfinity  1  Q1  3  O'clock 

6 

Inlet  Ptinfinity  2  Q1  9  O'clock 

7 

Pnozzle  static 

8 

9 

10 

11 

Kiel  Probe 

12 

Kiel  Probe 

13 

Kiel  Probe 

14 

Kiel  Probe 

15 

Kiel  Probe 

16 

Kiel  Probe 

17 

Kiel  Probe 

18 

Kiel  Probe 

19 

Kiel  Probe 

20 

Kiel  Probe 

21 

Kiel/TC  Combination 

22 

Kiel/TC  Combination 

23 

Kiel/TC  Combination 

24 

Kiel  Probe 

25 

Kiel  Probe 

26 

Kiel  Probe 

27 

Kiel  Probe 

28 

Kiel  Probe 

29 

Kiel  Probe 

30 

Kiel  Probe 

31 

Hub  Static  Pressures  P2  1 

32 

Hub  Static  Pressures  P2  2 

33 

Hub  Static  Pressures  P2  3 

34 

Hub  Static  Pressures  P2  4 

35 

Hub  Static  Pressures  P3  1 

36 

Hub  Static  Pressures  P3  2 

37 

Hub  Static  Pressures  P3  3 

38 

Hub  Static  Pressures  P3  4 

39 

Shroud  Static  Pressures  PI  3  O'clock 

40 

Shroud  Static  Pressures  PI  9  O'clock 

41 

Shroud  Static  Pressures  P3  3  O'clock 

42 

Shroud  Static  Pressures  P3  6  O'clock 

43 

Shroud  Static  Pressures  P3  9  O'clock 

44 

Shroud  Static  Pressures  Old  style  1  P3 

45 

Shroud  Static  Pressures  Old  style  2  P2 

46 

Shroud  Static  Pressures  Old  style  3  PI 

47i 

Shroud  Static  Pressures  Old  style  4  PI 

48 

Shroud  Static  Pressures  Old  style  5  Pi 

Time  References  are  from  turbine  out 

Table  3 .  Scanivalve  Port  Assignments . 
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The  logic  sequence  for  acquiring  data  from  the  48  ports 
of  the  Scanivalve  was  as  follows: 

1.  Reset  to  port  1 

2 .  Ensure  the  Scanivalve  was  actually  at  port  1 

3 .  Take  the  transducer  reading  using  the  DVM 

4.  Step  to  the  next  port 

5.  Ensure  the  Scanivalve  was  actually  at  the  next  port 

6.  Take  the  transducer  reading  using  the  DVM 

7 .  Go  to  step  4  until  all  48  ports  were  read 

2.  Sceuming  Digital  Voltmeter 

Connecting  a  16-channel  multiplexer  to  the  DVM  using 
the  'Analog  Bus'  allowed  the  multiplexer  and  DVM  to  operate 
together  as  a  Scanning  Multimeter.  ['Scanning  Multimeter' 
in^jlies  that  one  or  more  multiplexers  were  used  with  the 
multimeter  as  part  of  the  same  instrument  using  the  same  HP- 
IB  secondary  address]  .  The  power  of  this  configuration  was 
that  the  multimeter  was  extended  from  having  one  set  of 
terminals  to  having  sixteen  for  each  connected  multiplexer. 
Two  types  of  scanning  multiplexer  were  employed.  The  first 
type  was  a  thermocouple  relay  multiplexer  module  (HP  1347A) . 
A  thermistor  was  contained  inside  the  terminal  module,  and 
its  resistance  used  to  determine  temperature  with  which  to 
compensate  for  temperature  measurements  made  with  the 
thermocouples.  The  second  type  was  a  standard  relay 
multiplexer  module  (HP  1345A) ,  which  allowed  the  measurement 
of  sixteen  different  channels  using  a  single  multiplexer/DVM 
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connection.  The  channel  assignment  for  the  two  multiplexers 
is  given  in  Table  4 . 


Multiplexer 

Channel 

Thermocouples /Mux  100 

DVM  Scanner /Mux  200 

00 

Scanivalve  1 

Transducer 

01 

T^  Station  Inf  @  3 
o'clock 

Scanivalve  2 

Transducer 

02 

T^  Station  Inf  @  9 
o'clock 

Scanivalve  3 

Transducer 

03 

T^  Station  3  Probe  21 

Scanivalve  4 

Transducer 

04 

T^  Station  3  Probe  22 

05 

T,  Station  3  Probe  23 

06 

07 

08 

09 

10 

AT^  -  Station  Inf  @  3 
o'clock  and  Probe  21 

11 

AT^  -  Station  Inf  @  3 
o'clock  and  Probe  23 

12 

13 

14 

15 

Table  4.  Scanning  Multiplexer  Channel  Assignments. 


Readings  from  the  four  controlled  Scanivalve 
transducers  and  the  torque- flexure  strain  gauges  were  read 
directly  and  controlled  by  the  DVM  driver  software. 

To  measure  tenperature  directly,  the  thermocouple  wires 
were  connected  directly  to  the  thermocouple  terminal  module 
and  the  voltage  measured  and  converted  to  degrees 
Fahrenheit.  To  measure  the  instantaneous  differential 
tenperature  between  stations  3  and  Inf,  a  different  process 
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was  required.  The  thermocouples  were  connected  as  shown  in 
Fig.  8. 

The  logic  process  was  as  follows: 

1.  Measure  the  temperature  at  station  Inf  (using  the 
thermocouple  multiplexer) . 

2.  Convert  the  temperature  measured  in  step  1  to  volts 
using  the  equation  relating  temperature  to  voltage  for  a  J 
type  thermocouple. 

3.  Measured  the  voltage  between  the  two  thermocouples 
(using  the  scanning  multiplexer) . 

4.  Added  the  measured  voltage  of  step  3  to  the  voltage 
obtained  in  step  2 . 

5.  Converted  the  voltage  obtained  in  step  2  to  a 
tenperature  using  the  equation  relating  voltage  to 
tenperature  for  a  J  type  thermocouple. 
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3. 


Counter/Totali zer 


The  Counteir/ Totalizer  was  used,  to  measure  the  freguency 
produced  by  the  magnetic  speed  pickup.  The  pickup  produced 
six  pulses  per  revolution.  Therefore,  the  frequency 
measured,  in  Hz,  multiplied  by  10  produced  the  RPM  of  the 
rotor.  The  signal  produced  was  extremely  noisy  near  the  zero 
volt  level.  To  reject  this  noise,  the  level  to  trigger  the 
counter  was  set  to  +1  volt.  The  rotational  velocity  of  the 
rotor  was  verified  using  a  calibrated  strobe  light. 

C.  SOFTWARE 

HPVEE  was  the  software  used  to  control  the 
instrumentation.  Software  drivers  were  available  in  HPVEE 
for  all  instrximents  inside  the  mainframe.  The  'panel'  driver 
for  the  specific  instrument  served  two  purposes.  First,  it 
allowed  the  definition  of  a  measurement  state,  that  is  the 
settings  of  all  instrument  functions.  When  the  panel  driver 
operated,  the  corresponding  physical  instrument  was 
automatically  programmed  to  match  the  settings  defined  in 
the  panel  driver.  Secondly,  the  panel  driver  acted  as  a 
virtual  control  panel  for  interactively  controlling  the 
instrument. 

In  addition  to  controlling  the  instruments,  HPVEE  also 
allowed  immediate  reduction  of  the  measured  values  to 
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engineering  units.  This  allowed  immediate  feedback  of  data 
and  made  troubleshooting  simpler. 


D.  DATA  REDUCTION 

The  data  were  reduced  using  HPVEE.  In  the  description 
which  follows,  the  final  form  of  the  data  reduction  program 
is  described.  The  results  presented  were  calculated  with  the 
following  methods. 

The  mass  flow  rate  (in  Ibm/sec)  through  the  flow 
nozzle,  was  given  by 


m  =  3.3135 


C-  I — 

V  459.7 +  r, 


I  inf 


1-0.021926  -^ 


(1) 


where  /^)  was  in  inches  of  water  and  the  was  in 

inches  of  mercury.  The  constant  C  represented  the  discharge 
coefficient  and  was  calibrated  to  be  1.03. 

The  first  measure  of  the  power  absorbed  by  the 
compressor  was  obtained  from  the  measurements  of  torque,  M, 
and  speed,  N,  as 

HPci  =15865x10-^  M  N  (2). 

The  second  measure  of  power  was  obtained  from  the  flow 
rate  and  ten^erature  rise  across  the  stage,  namely 

HPc2  =  14137  Cp  m  lr, 3  -T,^)  (3)  . 

The  performance  of  the  coirpressor  was  described  in 
teirms  of  referred  quantities  that  retain  their  original 
\inits : 
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.  .  s 

mref=fn  — 


N, 


HPr^ 


HP 

5S 


M 

s 


(4) 


T  P 

where  6  =  -j^,  and  8  =  -r^,  and  where  T^^j=518.7°R  and 

^ REF  PrEF 

P^^f=29.92  inches  of  mercury. 

The  performance  to  be  evaluated  was  of  the  rotor  and 
stator  combination.  The  honeycomb  flow  straightener  in  the 
outlet  duct  was  necessary  in  order  for  the  torque  to  be 
measured,  but  the  losses  occurring  in  this  section  are  not 
included  in  appraising  the  conpressor  performance.  The 
total-to-total  (T-T)  pressure  ratio  and  efficiency  from 
station  1  to  station  3  were  the  parameters  which  gave  the 
performance  of  the  stage. 

The  total-to-total  efficiency  was  defined  as  the  ratio 
of  the  power  required  in  an  isentropic  conpression  to  the 
same  outlet  pressure,  to  the  actual  power  required.  The 
actual  power  was  given  by  equation  (2),  using  measurements 
of  torque  and  speed,  or  by  equation  (3),  using  measurements 
of  mass  flow  rate  and  terrperature  rise.  The  ideal  power  was 
given  by 
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mp,  =ui37  m.c,  ((r.^ (s) 

where  is  related  isentropically  to  the  pressures  by 


ftell 

r  p  ^ 

T 

imf  j 

J 

r-i 

r 


(6) 


where  y  is  the  ratio  of  specific  heats.  Substitution  into 
equation  (5) ,  gives 


f  T>  \ 


r-1 


43 


\  J 


-1 


(7)  . 


Three  different  methods  were  used  to  calculate  the 
total-to-total  efficiency.  Method  1  siitply  used  the 
temperature  rise  and  total-to-total  pressure  ratio  across 
the  stage 


IHPc 

HPc2 


^3  ^inf 


(8) 


The  second  method  used  inlet  static  pressure 
measurements  and  mass  flow  rate  to  calculate  an  inlet 
stagnation  pressure .  Defining 


where  the  area  corresponded  to  where  the  measurement  was 
taken,  it  follows  that 
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p  ^ 

r 

^>1 

\  ^inf  / 
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r-i 


’  2  * 


(10) 


where  and  V,;^  =  ^2  •  •  r,i„f  .  was  given  by 

^/inf 


X„=VV+1-A 

Thus 


inf 


(11)  . 
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inf 
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(12) 


was  an  efficiency  which  involved  the  mass  flow  rate 
measurement,  whereas  T|i  did  not. 

The  final  method  for  calculating  efficiency  depended  on 
the  measurement  of  torque  and  mass  flow  rate  but  did  not 
depend  on  the  measurement  of  the  temperature  rise  across  the 
stage,  namely 

(13)  . 


IHPc 

^=777^  =  ^- 
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ZV.  RESULTS  OF  THE  TEST  PROGRAM 


A.  INTRODUCTION 

The  perfoimance  of  the  stage  was  measured  at  50%,  60%, 
65%,  70%,  and  80%  of  the  design  speed  of  27,085  rpm.  The 
instrumentation  remained  constant  during  all  runs  and  no 
changes  to  the  data  acquisition  or  reduction  system  were 
required.  Some  changes  were  made  in  how  the  reduced  data 
were  saved  making  plots  of  torque,  horsepower  and  stage 
terrperature  rise  versus  mass  flow  rate  only  available  for 
70%  speed  and  80%  speed.  Each  set  of  plotted  experimental 
performance  data  was  fitted  with  a  second  order  least- 
squares  curve.  The  exit  impact  pressure  profiles  and  axial 
static  pressure  distributions  were  also  plotted. 

B.  STAGE  PERFORMANCE  MAPS 

The  total-to-total  pressure  ratio,  efficiency 
calculated  by  equations  (8),  (12)  and  (13),  referred  torque, 
referred  horsepower  and  stage  temperature  rise  were  plotted 
versus  the  referred  mass  flow  rate,  for  constant  speed 
conditions  expressed  as  a  percentage  of  design  speed. 

Figure  9  shows  the  total-to-total  pressure  ratio  vs 
referred  mass  flow  rate.  The  plot  shows  the  characteristic 
increase  of  pressure  ratio  and  decrease  of  mass  flow  rate  as 
the  flow  was  throttled.  At  all  speeds  except  for  80%  speed 
the  stage  was  throttled  to  stall.  The  loss  of  the  stage 
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occurred  before  stall  was  reached  at  80%  speed.  The  open- 
throttle  'choke'  condition  for  the  stage  could  not  be 
reached  due  to  the  pressure  drop  across  the  honeycomb  flow 
straightener . 

Figures  10,  11,  and  12  are  efficiency  plots.  The 
efficiency  T|,  was  expected  to  be  the  most  accurate  as  it  used 
torque  and  speed  measurements  which  provided  an  integral 
measurement  of  the  work  to  the  compressor.  The  other  two 
efficiency  measurements  used  the  tenperature  rise  across  the 
stage.  Since  the  total  temperature  was  not  expected  to  be 
constant  radially  at  the  exit  and  the  measurement  was  taken 
at  only  three  points  in  the  center  of  the  passage,  the  T|j 

and  Tij  measurements  of  efficiency  were  likely  to  be 
optimistic.  The  results  show  that  Tjj  varied  in  the  same  way 
as  the  other  two  efficiency  measurements,  but  was 
consistently  lower  in  magnitude.  The  data  for  rjj  at  80% 

speed  measured  2.5%  lower  than  the  other  constant  speed 
lines.  This  was  possibly  the  result  of  drift  in  the  strain- 
gauge  bridge  output  (caused  by  increased  tenperature )  which 
was  not  balanced  to  zero  before  the  final  run  was  started. 
The  drift  was  observed  after  the  loss  of  the  stage.  The 
results  forT|j  and  support  this  contention  since  no 
appreciable  drop  was  observed  in  the  80%  speed  results. 
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The  CFD  analysis  results  shown  on  the  performance  maps 
were  obtained  for  the  rotor  geometry,  and  did  not  include 
the  stator.  These  results  are  discussed  in  the  following 
section. 

Figure  13  shows  the  increased  torque  required  as  the 
coitpressor  was  throttled.  Figure  14  shows  the  increase  in 
the  stage  tenperature  rise  as  the  coitpressor  was  throttled. 
The  stage  tenperature  rise  was  calculated  as  the  average  of 
the  differential  tenperature  of  one  probe  at  station  Inf 
with  two  probes  at  station  3 .  Figure  15  shows  a  comparison 
of  the  horsepower  required  to  drive  the  stage  calculated 
using  the  mass  flow  rate  and  stage  temperature  rise,  (Eq. 
(3)),  with  that  calculated  using  the  torque  and  speed 
measurements  (Eq.  (2)).  The  horsepower  measured  by  torque 
and  speed  was  consistently  greater  than  horsepower  deduced 
from  mass  flow  rate  and  temperature  rise. 

C.  PRESSURE  PROFILES 

Figures  16  and  17  show  the  exit  total  pressure  profiles 
at  70%  speed  and  80%  speed  respectively.  (Note  that  the 
probes  were  at  different  peripheral  settings  in  the  radial 
direction  shown  in  Table  2)  .  The  ratio  of  each  measured 
total  pressure  at  station  3  to  the  total  pressure  at  station 
1  was  calculated.  Three  different  throttle  settings  are 
shown  at  each  speed.  At  design  conditions,  the  stage  was 
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designed  to  have  nearly  uniform  total  pressure  across  the 
exit  [Ref.  1].  The  plots  show  that  as  the  stage  was 
throttled,  at  an  intermediate  throttle  setting  the  total 
pressure  was  nearly  uniform  across  the  annulus,  except  near 
the  walls.  The  plots  also  show  that  as  stall  was  approached, 
which  corresponded  to  the  lowest  mass  flow  rate  curve,  a 
separation  had  progressed  radially  outward  from  the  hiib. 

Figures  18  and  19  are  plots  of  the  static  pressure 
measured  at  different  axial  stations  at  the  same  rotor 
speeds  and  throttle  settings  shown  in  the  data  in  Figs.  16 
and  17  respectively.  All  pressures  are  scaled  by  the 
pressure  at  station  infinity. 

At  station  2  a  strong  static  pressure  gradient  existed 
from  hub-to-tip  at  all  measured  states  of  throttling  and 
wheel  speed.  This  gradient  is  consistent  with  the  high 
tangential  velocities  expected  at  the  rotor  exit.  The 
pressure  rise  through  the  stator  behaved  in  a  consistent  way 
with  throttling  at  each  of  the  speeds  measured.  First,  the 
pi^sssure  rise  at  the  hub  was  large  at  all  throttle  settings . 
At  the  tip,  the  rise  improved  with  closing  throttle,  xintil 
the  static  pressure  difference  between  hub  and  tip  at 
station  3  disappeared  as  the  stage  approached  stall. 
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Pressure  Ratio  (T-T)  vs  Referred  Mass  Flow  Rate 


Adiabatic  Efficiency 


Adiabatic  Efficiency 


Referred  Mass  Flow  Rate  (Ibm/sec) 


Fig.  12  Adiabatic  Efficiency  (t|,)  vs  Referred  Mass  Flow 
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Fig.  13  Referred  Torque  vs  Referred  Mass  Flow  Rate. 


Fig.  15  Referred  Horsepower  vs  Referred  Mass  Flow  Rate. 


41 


Radial  Distance  y  (inches) 


Hub 

2r 


Fig.  16  Exit  Pressure  Distribution  at  70%  Speed. 
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Radial  Distance  y  (inches) 


Hub 


Stage  Total  Pressure  Ratio 


Fig.  17  Exit  Pressure  Distribution  at  80%  Speed. 


43 


Fig.  18  Axial  Static  Pressure  Distribution  at  70%  Speed, 
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Pstatic/PInf 


Fig.  19  Axial  Static  Pressure  Distribution  at  80%  Speed. 
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V.  RESULTS  OF  CFD  ANALYSIS 


A.  INTRODUCTION 

Before  starting  the  experimental  work,  the  flow  through 
the  rotor  was  analyzed  and  the  rotor  performance  was 
predicted  using  CFD.  Specifically,  Rotor  Viscous  Code  3-D 
(RVC3D) ,  a  computer  code  for  analysis  of  three-dimensional 
viscous  flows  in  turbomachinery  [Ref.  3],  was  used  to 
compute  performance  maps  of  the  rotor  from  50%  to  100%  of 
the  design  speed.  Additionally,  using  the  graphics  package 
FAST  from  NASA  Ames,  a  number  of  parameters  were  calculated 
and  graphically  displayed  for  100%  speed.  The  predicted 
results  for  the  rotor  were  coitpared  to  the  experimentally- 
obtained  maps  for  the  stage  and  they  were  observed  to  be 
numerically  consistent  with  each  other. 

B.  (HlID  CSNERATION 

The  geometry  data  for  the  hub,  tip  and  blade  for  the 
stage  were  provided  electronically  by  NASA  Lewis.  The  data 
were  in  MERIDL  (z,  r,  theta,  dtheta)  format. 

The  geometry  file  was  entered  into  Turbomachinery  C 
GRID  (TCGRID)  which  generated  a  three-dimensional  C  type 
grid  cortpatible  with  the  Rotor  Viscous  Code  3-D  (RVC3D} 
[Ref.  7].  The  grid  which  was  185X40X49  (in  i,j,  and  k)  is 
shown  in  Fig.  20.  The  output  grid  was  stored  in  plot3d 
multiple  grid  format.  A  FORTRAN  program  was  used  to  rotate 
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the  grid  one  grid  width  to  form  the  multi-block  grid  that 
was  an  input  to  FAST. 


Fig.  20  Con^utational  Grid  for  RVC3D  Calculations 
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C.  COMPUTATION 

The  code  solved  the  thin-layer  Navier-Stokes  equations 
with  an  explicit  four  stage  Runge-Kutta  scheme.  The  number 
of  stages  was  optional;  however,  four  was  recommended  by  the 
author.  A  spatially-varying  time  step  (At)  was  used  to 
accelerate  the  convergence  of  the  code.  At  was  updated 
every  fifty  iterations. 

Since  the  flow  was  supersonic,  a  second-difference 
artificial  viscosity  term  was  used.  The  value  was  set  to  1.0 
in  order  to  produce  a  well  defined  shock  with  little 
smearing  or  waviness.  No  first-order  artificial  viscosity 
term  was  required. 

RVC3D  used  an  implicit  residual  smoothing  routine  which 
allows  the  multistage  scheme  to  run  at  a  high,  unstable 
Courant  number,  but  stability  was  maintained  by  smoothing 
the  residual  using  an  implicit  filter.  The  initial  Courant 
number  chosen  was  5.0. 

The  code  ran  conp»letely  non-dimensionally,  with 
variables  referenced  to  speed  of  soiand,  reference 
temperature  and  reference  pressure.  The  input  parameters 
that  controlled  the  flow  were  the  normalized  blade  row 
rotational  speed  per  unit  length  and  the  ratio  of  the  hiob 
exit  static  pressure  to  the  reference  total  pressure. 
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D.  RESULTS 

Version  208  of  the  code  was  run  on  three  different 
machines,  a  Silicon  Graphics  Indigo2  workstation,  the  NPS 
Cray,  and  NASA  Ames  Cray.  On  the  NASA  Ames  Cray,  each 
iteration  took  2 . 3  seconds  of  CPU  time  and  required  9Mw  of 
storage.  The  code  was  cortpiled  using  FORTRAN  90  and  02 
optimization  that  gave  good  scalar  and  vector  optimization 
and  performed  autotasking.  Six  processors  were  normally 
used.  In  order  for  the  grid,  solution,  and  restart  files  to 
be  compatible  with  each  machine,  IEEE  format  was  selected 
for  the  data  input  and  output. 

To  compute  the  operating  lines  at  a  constant  rotational 
speed,  from  the  choke  region  to  the  stall  region,  the  ratio 
of  the  hub  exit  static  pressure  to  the  reference  total 
pressure  {p^^J  was  increased  in  steps.  For  each  p^^,  at  which 
a  solution  converged,  the  mass  flow  rate,  the  mass-averaged 
and  exit  total  pressures,  and  the  adiabatic  efficiency 
were  recorded  from  the  output  file  for  plotting.  When  a 
converged  solution  could  not  be  obtained  (near  stall) ,  the 
rotational  speed  was  changed  and  a  new  range  of  values  of 
Prat  was  selected  and  entered.  The  operating  speeds  computed 
in  this  manner  were  70%,  80%,  90%  and  100%  of  design  speed. 

tionally ,  one  point  was  computed  at  50%  of  design  speed. 
The  results  are  shown  plotted  in  Figs.  21  and  22.  The 
results  were  observed  to  be  extremely  close  to  the  results 
Sanger  obtained  using  an  Euler  solver  with  boundary  layer 
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modeling  (Denton's  Code  [Ref.  1]).  At  design  mass  flow  rate 
of  17.1  Ibm/sec  Sanger  obtained  a  total  pressure  ratio  of 
1.61,  and  using  RVC3D  a  total  pressure  ratio  of  1.601  was 
found. 

The  predicted  results  for  the  rotor  at  70%  and  80% 
speed  are  shown  plotted  with  the  experimental  data  for  the 
stage  in  Figs  9  and  10.  A  direct  numerical  conparison  is  not 
meaningful  since  the  predicted  results  were  only  for  the 
rotor.  However,  the  conputed  results  are  qualitatively 
consistent  with  what  was  measured  for  the  stage  behavior. 
The  pressure  ratio  conputed  for  the  rotor  was  consistently 
higher  than  was  measured  for  the  stage,  and  the  range  of 
mass  flow  rate  predicted  at  each  rotor  speed  was  consistent 
with  the  stage  measurements .  At  80%  speed,  solutions  were 
obtained  closer  to  choke  (more  open  throttle)  than  could  be 
achieved  in  the  test  program  for  the  stage. 
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Rotor  Total  Pressure  Ratio 


Fig.  21  Rotor  Total  Pressure  Ratio  vs  Corrected  (Referred) 
Mass  Flow  Rate  Coinputed  Using  RVC3D 
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Fig.  22  Rotor  Adiabatic  Efficiency  vs  Corrected  (Referred) 
Mass  Flow  Rate  Confuted  Using  RVC3D 
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VI.  CONCLUSIONS  AND  RECOMHENDATIONS 


A  test  program  was  conducted  to  80%  of  design  speed  to 
evaluate  the  performance  of  a  new  transonic  axial  compressor 
stage.  Tests  to  100%  speed  were  prevented  by  the  failure  of 
the  spinner  retaining  bolt  which  led  to  the  loss  of  the 
stage.  The  cause  of  the  failure  was  shown  by  analysis  and 
tests  will  resxame  when  replacement  parts  are  delivered. 

A  new  data  acquisition  system,  using  PC-based  control 
of  VXI-bus  instrxamentation,  was  programmed  and  used 
successfully  in  the  tests.  HPVEE  software  was  used  with  the 
VXI-bus  hardware  for  Scanivalve  control,  acquisition  of 
pressures,  tenperatures ,  strains  and  frequencies,  and  for 
the  reduction  of  data.  The  system  was  foiind  to  be  extremely 
flexible  and  well-suited  to  the  researcher  in  the  laboratory 
environment . 

The  flow  through  the  rotor  alone  was  cortputed  using  the 
RVC3D  code.  The  test  results  obtained  agreed  qualitatively 
with  the  predicted  rotor  performance  and  were  also 
numerically  consistent  with  them.  This  was  a  first  step  in 
the  code  validation  process .  The  next  step  will  be  to 
measure  the  flow  profile  at  the  exit  of  the  rotor  and 
compare  what  is  measured  to  the  predicted  flow  field.  The 
final  step  will  be  to  measure  the  flow  between  the  rotor 
blades  and  pressure  distribution  on  the  blade  surfaces  and 
compare  these  to  predicted  results. 
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with  the  essperience  obtained  from  the  present  study  the 

following  specific  recommendations  are  made: 

•  With  the  inlet  temperature  fluctuations  caused  by  the 
coastal  location  of  the  laboratory,  increased  use  of 
averaging  is  required.  First,  a  scan  of  tenperatures  and 
wheel  speed  should  be  conducted,  then  a  pressure  scan, 
then  a  second  scan  of  tenperatures  and  wheel  speed  should 
be  taken  and  averaged  with  the  first  one  to  define  one 
data  point .  Three  data  points  should  then  be  taken  and 
averaged  to  produce  an  averaged  data  point  to  plot  on  the 
performance  maps. 

•  The  torque- flexure  strain  gauge  bridge  output  should  be 
readjusted  to  zero  between  each  constant  speed  nan.  (The 
output  should  be  zero  with  no  torque  applied.  However, 
the  zero  tended  to  drift  with  increased  operating 
tenperatures  and  therefore  needed  to  be  adjusted) . 

•  Install  a  strain  gauge  multiplexer  for  strain 
measurements.  Software  procedures  would  then  make  the 
previous  recommendation  aannecessaary  and  would  be  the 
preferred  solution. 

•  In  rebuilding  the  stage,  use  only  a  solid  steel  spinner 
retaining  bolt,  use  a  transparent  case  wall,  and  re-use 
the  video-monitoring  and  recording  system  demonstrated  in 
the  present  study. 


56 


LIST  OF  REFERENCES 


1.  Sanger,  N.  L. ,  Design  of  a  Low  Aspect  Ration  Transonic 
Compressor  Stage  Using  CFD  Techniques.  ASME  Paper  No.  94- 
GT-236,  Presented  at  the  International  Gas  Turbine  and 
Aeroengine  Congress  and  Exposition,  The  Hague,  The 
Netherlands,  June  13-16,  1994. 

2.  Reid,  William  D. ,  Transonic  Axial  Compressor  Design  Case 
Study  and  Preparations  for  Testing.  Master's  Thesis, 

Naval  Postgraduate  School,  Monterey,  California, 

September  1995. 

3.  Chima,  Roderick  V. , "RVC3D  -  Rotor  Viscous  Code  3-D  User's 
Manual  and  Documentation",  NASA  Lewis  Research  Center, 
March,  1992. 

4.  United  Sensor  Corporation,  "List  Price  Schedules",  3 
Northern  Boulevard,  Amherst,  NH  03031-2329,  April  04, 
1994. 

5.  Shreeve ,  R .  P .  ,  Report  on  the  Testing  of  a  Hybrid 
(Radial-To-Axial)  Compressor.  Technical  Report  NPS- 
57Sf73112A,  Naval  Postgraduate  School,  Monterey, 
California,  Ndvember  1973. 

6.  Geopfarth,  Robert  N. ,  Development  of  a  Device  for  the 
Incorporation  of  Multiple  Scanivalves  into  a  Computer- 
Controlled  Data  System.  Master's  Thesis  NPS67-79-002 , ■ 
Naval  Postgraduate  School,  Monterey,  California,  March 
1979. 

7.  Chima,  Roderick  V. ,  "TCGRID  -  Turbomachinery  C-grid 
User's  Manual",  NASA  Lewis  Research  Center  ,  November 
1990. 


57 


58 


APPENDIX  A  ANALYSIS  OF  THE  SPINNER  BOLT  FAILURE 

A.  OVERVIEW 

While  operating  at  80%  speed,  with  the  inlet  throttled 
to  produce  near  stall  conditions,  the  bolt  attaching  the 
spinner  to  the  rotor  failed.  The  spinner  detached  and 
contacted  the  rotor,  causing  complete  loss  of  the  stage  and 
casing.  Photos  of  the  failed  bolt,  the  destroyed  stage,  and 
the  damaged  Lucite  casing  are  shown  in  Figs.  Al,  A2,  and  A3 
respectively.  While  running  the  coitpressor,  a  video  camera 
was  being  used  to  monitor  and  film  the  operation.  Using  the 
video  film  (at  30  frames/ sec)  and  pressure  readings  measured 
immediately  before  the  event,  an  analysis  was  conducted  of 
the  failure. 

B.  CONSTRUCTION 

The  spinner  used  with  the  stage  had  been  modified  from 
the  design  used  on  the  Vavra  stage.  The  retaining  bolt  was 
required  to  be  two  inches  longer;  namely,  3.09  inches  versus 
1.09  inches  for  the  Vavra  spinner  retaining  bolt.  The  bolt 
was  a  Socket  Head  Cap  Screw  #8-32,  which  was  not 
manufactured  (standard)  with  the  new  required  length.  In 
order  to  increase  the  length  of  the  bolt,  a  1-1/8  inches 
bolt  was  cut  in  half.  A  metal  sleeve  with  5/32  inches  o.d. 
and  1/8  inches  i.d.  was  placed  over  the  two  halves  that  had 
been  turned  down  to  1/8  inches  o.d.  The  metal  sleeve  was 
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then  silver  brazed  at  1050°  F  to  form  the  extension.  After 
spin  balancing  the  rotor  and  drive  shaft,  the  spinner  was 
attached  using  the  bolt  and  the  assembly  rebalanced.  This 
required  removing  material  from  the  inner  edge  of  the 
spinner.  When  balance  was  finally  achieved,  the  bolt 
attaching  the  spinner  was  'lightly  torqued' .  Only  after  the 
event  was  it  found  that  the  metal  was  cut  from  copper 
tubing . 

C.  ANALYSIS 

In  order  to  estimate  the  aerodynamic  load  on  the 
spinner  the  pressure  measurements  recorded  just  before  the 
failure  occurred  were  used.  The  throttle  had  been  closed  a 
small  fraction  more  when  the  failure  occurred;  therefore, 
the  previously  scanned  pressure  values  are  considered  to  be 
very  close  to  those  that  were  present.  For  a  worst-case 
analysis,  the  pressure  inside  the  spinner  was  assumed  to  be 
the  stagnation  pressure  at  the  nose  of  the  spinner.  (This 
asstimption  is  r:easonable  since  the  hole  that  the  bolt  passed 
through  was  not  sealed  and  leakage  in  would  have  occurred. 
The  shoulder  at  the  rim  of  the  spinner  was  also  not  sealed 
and  the  pressure  behind  the  spinner  would  leak  out  and  be 
less  than  this  'worst  case'  assumption) .  The  pressure  on  the 
outer  wall  was  taken  to  be  the  average  of  the  static 
piressures  taken  at  station  one,  another  static  pressure 
reading  taken  mid-length  of  the  spinner,  and  the  static 
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pressure  at  the  nose  of  the  spinner.  All  the  static  pressure 
measurements  were  taken  at  the  casing  wall.  With  this 
pressure  distribution  an  axial  force  of  11.0  Ibf  was 
calculated.  This  axial  force  would  have  caused  a  tensile 
stress  of  1.6  ksi  in  the  copper  sleeve  of  the  bolt-  The 
stress  in  the  copper  sleeve  due  to  the  torque  used  to  attach 
the  bolt  to  the  rotor  was  estimated  to  have  been  10.3  ksi. 
The  maximum  steady  total  tensile  stress  in  the  copper  sleeve 
was  therefore  11 . 9  ksi . 

From  the  videotape  of  the  event,  it  was  observed  that 
the  spinner  went  upstream  and  disappeared  from  view  within  2 
video  frames,  came  back  and  'flew'  stably  for  several 
seconds,  moving  arovind  the  annulus,  and  then  fell  into  the 
blading.  Since  the  frame  speed  was  known  to  have  been  30 
frames /sec,  an  estimate  of  the  acceleration  of  the  spinner 
as  it  left  the  rotor  was  made  by  comparing  the  relative 
positions  of  the  spinner  over  two  frames.  The  estimated 
average  acceleration  was  225  ft/sec"".  The  drag  force  on  the 
spinner  (downstream)  was  calculated  to  have  been  2.52  Ibf. 
The  net  force  required  to  accelerate  the  spinner  upstream 
was  calculated  to  have  been  2.33  Ibf.  The  initial  force  on 
the  spinner  necessary  to  explain  the  video  images  was 
therefore  4.85  Ibf.  This  coitpares  to  the  estimate  made  using 
pressure  measurements  that  calculated  a  force  of  11.0  Ibf . 
The  image  on  the  first  frame  following  detachment  was 
blurred  as  the  spinner  moved  rapidly  upstream,  making  an 
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accurate  estimation  of  the  instantaneous  acceleration  at  the 
time  of  detachment  impossible.  Therefore,  although  the 
calculations  do  not  agree  closely,  they  do  demonstrate  that 
the  forces  calculated  to  be  present  were  of  the  order  of 
magnitude  necessary  to  explain  the  observed  behavior  of  the 
spinner  if  the  copper  sleeve  failed  in  tension. 

The  tensile  strength  of  annealed  copper  is  32  ksi  and 
the  yield  strength  is  5  ksi.  Copper  has  the  property  of 
strain  hardening.  Therefore,  each  time  the  sleeve  was 
stressed  into  plastic  deformation  the  tensile  strength  would 
have  decreased.  With  a  stress  twice  the  yield  strength 
present,  failure  of  the  bolt  occurred  when  the  tensile  force 
exceeded  the  reduced  tensile  strength. 
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Fig.  Al  View  of  the  Failed  Spinner 
Retaining  Bolt 


Fig.  A3  View  of  Lucite  Casing  After  Failure 
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APPENDIX  B  RESULTS 


The  reduced  results  of  all  test  runs  are  presented 
below.  Data  points  with  no  values  present  for  horsepower 
were  caused  by  problems  with  the  data  acquisition  system. 
The  point  at  the  same  state  was  re-measured  until  a  complete 
data  scan  was  obtained.  The  corresponding  raw  data  file 
raw_tcr.dat  is  not  printed  but  is  held  on  a  3.5  inch  floppy 
disk  by  Dr.  Shreeve. 


67 


63005636295  7.74  9.01  95,72  108.02  90  65  10?  :^n  17M9R4  .ifl  ~~r7i  TTSa - KTK - - oHz 

63005636491  7.21  8.45  95.58  109,41  89.10  101,99  17474,72  394,64  1,22  0.85  '~0.'87  ~0li0  6?52 

63005636725  7,49  8.74  94,98  108.05  89.94  102.32  17418.48  390.98  1.22  0.87  ~a8e  612  643? 

630056369561  8.21  |  9.44  |  93.29  |  104.40  |  91.76  |  102.69  |  17367.271  378.89  |  1.21  |  0.89  0.90  ~0l8  6^ 


^?^f*^f^:^:^5^5^<c>5^50  50  50CD<o<DtDcoocDCDcDS<o^^<D<D®^ 
n->i<Doairooio>*^03^<D->^03-vjocoXro^ooia>^5S5S^2£S2.® 


a  (Q  3 
3 


Referrei 


Referred  HP  HP 

Data  Mass  Mass  HP  -  Referred  HP  -  Referred  -  RPM  Torque  Pressure 


APPENDIX  C  SOFTWARE  USER'S  GUIDE 


A.  HPVEE 

The  procedure  for  operating  the  HPVEE  data  acquisition 
routine  was  as  follows : 

1.  The  input  array  for  the  Scani valve  Control  user's 
object  was  initialized.  The  0"“  element  was  the 
Scanivalve  Position  (1  through  4)  and  the  next  48 
elements  of  the  array  indicated  whether  the 
corresponding  port  should  be  measured  (0  -  no  skip, 
1  -  yes  measure) 

2.  When  the  rig  was  stable  and  at  the  desired  condition 
for  measurement  the  START  button  was  depressed. 
Desired  quantities  were  input  as  prompted. 

3 .  When  the  data  acquisition  process  was  finished  the 
data  were  automatically  written  to  the  hard  disk. 

4.  The  operating  condition  was  changed  and  the  process 
repeated  from  step  1. 

All  programs  used  were  saved  on  a  3.5  inch  floppy  and 
are  in  the  possession  of  Dr.  Shreeve.  A  graphical  listing  of 
the  program,  with  redundant  portions  for  scanivalve 
positions  2  through  4  removed,  follows: 
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C:\BARTMcr_recluction_print.vee 


at:  C5  AlphaNumecic  "SV  in  volts” 

at;  D4  Pormula  "Pressure  ins  Hg  absolute 

Formula:  {A-A(  1 1 )  *  (8/ (All!) -A(l  1 )  )+C 
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Main:  Temperature  Cates:  Panel  View 

100  200  300  400 


PanelView  "Pressure  Measurement"  (w=1000,  h=802) 
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Main:  Pressure  Measurement:  Panel  View 
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Sat  20/Sep/1997  12:23:32 

Main;  Preaaure  Measurement;  Scanivalve  Control  Object  (M.197,4);  Detail  View,  Total  Objects:  10 
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<42  AND  N>3 


C:\BART\tcr_reduction_print.vee 

Sat  20/Sep/1997  12:23:32 

Main:  Pressure  Measurement:  Scanlvalve  Control  Object:  SV 1  Data  Acquisition:  CHECK  THE  PORT  #  SV1:  CONVERT^TO^BCD  (M.197.4.45.1.2):  Detail  View.  Total  Objects:  7 
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C:\BART\tcr_reduction_prlnt.vee 

Sat  20/Sep/1997  12:23:32 

Main:  Pressure  Measurement:  Scanlvalve  Control  Object:  SV 1  Data  Acquisition:  Step  SV1  (M.197.4.45.13):  Detail  View.  Total  Objects:  7 
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Formula:  {A+B)/2 
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CABARTMcr^reductionjjrint.vee 
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AlphaNumeric  "pnozzle"  (x=150,  y«10)  (w=135,  h=50) 
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C:\BARTMcr_reduction_print.vee 


B.  RVC3D 

All  software  used  and  results  obtained  using  RVC3D  and 
associated  programs  are  stored  in  a  directory 
/ jedi/d3/blgrossm  on  the  NFS  Cray.  Dr.  G.  V.  Hobson  has 
access  to  the  files.  The  procedure  for  operating  RVC3D  was 
as  follows: 

1.  Input  the  geometry  provided  as  geometry.dat  into 
TCGRID  to  obtain  the  grid  file  fort.l.  geometry. out  was 
produced  and  examined  to  ensure  the  grid  was  properly 
produced.  The  command  used  was  "tcgrid.exec  <geometry . dat 
>geometry . out" . 

2.  Modified  rvc.in  to  reflect  the  desired  initial 
conditions  to  be  analyzed  [Ref.  3] . 

3.  Submitted  the  job  to  the  queue  using  rvc3d.nqs. 

4.  After  conpletion,  a  solution  file  q.sol  and  a 
result  file  rvc.out  were  output. 

5.  If  further  iterations  were  desired,  modified  rvc.in 
such  that  iresti=l  and  renamed  q.sol  to  q.res. 

6.  Resvibmitted  the  job  to  the  queue. 

Anple  assistance  was  available  from  computer  center 
personnel.  Examples  of  geometry.dat,  rvc.in  and  rvc3d.nqs 
are  included  in  the  following  pages. 
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Listing  of  rvcSd.nqs: 


#QSUB  -q  reg  -IT  15280  -It  15280  -IM  9Mw  -me 
set  RUN  =  grid 

assign  -F  til  -N  ieee  -a  ${RUN}.ieee  u:l 
set  RUN  =  q 

assign  -F  til  -N  ieee  -a  ${RUN}.sol  u:3 
assign  -F  til  -N  ieee  -a  ${RUN}.res  u:2 


cd  /d3/blgrossm/tentproj 
rvc3d  <rvc.in  >rvc.out 


Listing  of  rvc.in: 

■sanger  rotor* 

toll  im=185  jm=40  ]an=49  itl=20  iil=85  &end 

tol2  nstg=4  cfl=5  aviscl=0  avisc2=1.0  avisc4=.5  irs=l  epi=0.60  epj=0.70 
epk=0.70  ivdt=l  itmax=5100  tend 
tol3  ibcin=l  ibcex=3  ires=10  icmt=50  iresti=l  tend 
tol4  om=-2.54  prat=1.1905  emxx=.46  alex=52  igeom=l  S^end 
tol5  ilt=2  iltin=2  renr=6.651e6  tw=0  vispwr=.6667  crauitm=14  jedge=25 
kedge=18  srtip=0  xrle=-.2  xrte=0.3  dblh=0.001  dblt=0.02  &end 
tol6  &end 


Listing  of  geometiry.dat : 

toaml  merid=3  im=185  jm=40  ]aa=49  itl=20  icap=16  k2d=ll  &end 
toain2  nle=15  nte=10  dsle=3.5e-4  dste=2.5e-4  dshub=1.5e-5  dstip=1.5e-5 
dswte=3.5e-4  dswex=.010  dsthr=.39  dsmin=1.2e-5  dsinax=.002 
dsin=.010  clear=. 00000  rcom=.008  &end 
toaiti3  itenn=100  idbg=0  0  0  0  0  0  0  0  aabb=.6  ccdd=.2  &end 
toam4  zbc=-0.100  -0.040  0.2000  -0.1000  -0.0400  0.2000 

rbc=  0.220  0.220  0.2800  0.4560  0.4560  0.4560  &end 

■NPGS  TRANSONIC  DESIGN  -GENER.  FROM  V8FBND13 INCDP-BASE=DEC13  INNER* 

32  14 

-.833  -.667  -.5833  -.50  -.4583  -.425  -.40833  -.3917 
-.3333  -.25  -.16666  -.0833  0.  .0833  .1  .1167 

.125  .133  .14167  .15  .1583  .16666  .175  .183 
.19166  .208  .25  .3333  .41666  .5  .667  .8333 

.00083  .00083  .00083  .00083  .001667  .00625  .01083  .01667 
.04467  .08937  .13405  .17873  .2234  .2681  .277  .286 
.2904  .29333  .295583  .297667  .29933  .30033  .30166  .3025 
.303333  .304166  .306083  .306083  .306083  .306083  .306083  .306083 

-.8333  -.7637  -.666666  -.5  -.33333  -.2083  -.1239  -.04167 
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.0529  .2167  .3083  .5  .66666  .83333 


.4583  .4583  .4583  .4583  .4583  .4583  .4583  .4583 
.4583  .4583  .4583  .4583  .4583  .4583 


11  23  22 

0.000000  0.000514 
0.025568 

0.031878  0.038395 
0.081719 

0.089506  0.097415 
0.001360  0.001751 
0.026281 

0.032350  0.038563 
0.079987 

0.087529  0,095222 
0.002013  0.002353 
0.026682 

0.032644  0.038731 
0.078992 

0.086358  0.093887 
0.002617  0.002913 
0.027000 

0.032848  0.038800 
0.077748 

0.084894  0.092212 
0.003875  0.004123 
0.027672 

0.033327  0.039064 
0.075885 

0.082643  0.089584 
0.007383  0.007577 
0.029584 

0.034782  0.040025 
0.072754 

0.078711  0.084862 
0.010925  0.011094 
0.031430 

0.036196  0.040990 
0.070365 

0.075575  0.081001 
0.012591  0.012749 
0.032161 

0.036707  0.041277 
0.069210 

0.074074  0.079161 
0.014092  0.014243 
0.032945 

0.037327  0.041734 
0.068699 

0.073344  0.078207 
0.014068  0.014213 
0.032800 

0.037163  0.041553 
0.068472 

0.073071  0.077861 


0.002301 

0.045135 

0.104086 

0.003532 

0.044956 

0.101740 

0.004135 

0.044945 

0.100278 

0.004692 

0,044857 

0,098434 

0.005878 

0.044881 

0.095503 

0.009239 

0.045312 

0.090132 

0.012640 

0.045811 

0.085687 

0.014226 

0.045871 

0.083566 

0.015665 

0.046166 

0.082428 

0.015624 

0.045971 

0.082024 


0.005027 

0.052088 

0.109470 

0.006238 

0.051567 

0.107019 

0.006837 

0.051351 

0.105463 

0.007384 

0.051036 

0.103490 

0.008529 

0.050778 

0.100323 

0,011743 

0.050643 

0.094442 

0.014965 

0.050659 

0.089541 

0.016447 

0.050490 

0.087197 

0.017803 

0.050623 

0.085913 

0.017747 

0.050*417 

0.085465 


0.008746 

0.059237 

0.113531 

0.009908 

0.058387 

0.111013 

0.010492 

0.057968 

0.109392 

0.011017 

0.057417 

0.107325 

0.012096 

0.056765 

0.103986 

0.015099 

0.056018 

0.097727 

0.018075 

0.055534 

0.092492 

0.019417 

0,055133 

0.089983 

0.020663 

0.055105 

0.088589 

0.020588 

0.054890 

0.088111 


0.013525 

0.066568 

0.116248 

0.014595 

0.065405 

0.113692 

0.015144 

0.064787 

0.112029 

0.015627 

0.064001 

0.109902 

0.016608 

0.062939 

0.106450 

0.019322 

0.061437 

0.099942 

0.021980 

0.060436 

0.094488 

0.023145 

0.059800 

0.091869 

0.024255 

0.059611 

0.090402 

0.024157 

0.059390 

0.089904 


0.019448 

0.074067 

0.116636 

0.020362 

0.072609 

0.114041 

0.020848 

0.071798 

0.112356 

0.021259 

0.070780 

0.110207 

0.022099 

0.069314 

0.106724 

0.024430 

0.066995 

0.100171 

0.026692 

0.065366 

0.094690 

0.027641 

0.064492 

0.092059 

0.028587 

0.064143 

0.090584 

0.028465 

0.063918 

0.090081 


100 


0.014718 

0.032972 

0.037239 

0.067997 

0.072518 

0.222890 

0.236541 

0.239910 

0.266519 

0.270677 

0.251549 

0.262257 

0.264864 

0.285332 

0.288573 

0.267527 

0.276842 

0.279093 

0.296593 

0.299375 

0.284586 

0.292532 

0.294438 

0.309070 

0.311399 

0.309498 

0.315614 

0.317067 

0.328005 

0.329743 

0.357485 

0.360641 

0.361380 

0.366779 

0.367626 

0.398682 

0.400090 

0.400417 

0.402763 

0.403120 

0.418757 

0.419585 

0.419778 

0.421154 

0.421360 

0.428363 

0.428980 

0.429124 

0.430151 

0.430303 

0.441651 

0.442033 

0.442122 

0.442760 

0.442854 


0.014854 

0.016226 

0.018291 

0.021057 

0.024535 

0.028737 

0.041537 

0.045868 

0.050230 

0.054624 

0.059050 

0.063508 

0.077153 

0.081192 

0.084538 

0.087114 

0.088863 

0.089032 

0.223164 

0.224119 

0.225574 

0.227560 

0.230111 

0.233273 

0.243389 

0.246988 

0.250699 

0.254516 

0.258430 

0.262434 

0.274899 

0.251717 

0.278461 

0.252483 

0.281335 

0.253645 

0.283503 

0.255222 

0.284954 

0.257236 

0.285162 

0,259714 

0.267534 

0.270281 

0.273121 

0.276052 

0,279067 

0.282162 

0.291878 

0.267655 

0.294678 

0.268328 

0.296946 

0.269348 

0.298663 

0,270728 

0.299814 

0.272485 

0.299964 

0.274639 

0.281391 

0.283738 

0.286156 

0.288655 

0.291230 

0.293877 

0.302218 

0.284683 

0.304631 

0.285263 

0.306589 

0.286140 

0.308072 

0.287324 

0.309068 

0.288826 

0.309192 

0.290661 

0.296378 

0.298351 

0.300365 

0,302444 

0.304590 

0.306799 

0.313784 

0.309561 

0.315811 

0.310012 

0.317459 

0.310694 

0.318709 

0.311611 

0.319549 

0.312770 

0.319648 

0.314181 

0.318542 

0.320037 

0.321552 

0.323091 

0.324678 

0.326317 

0.331527 

0.357513 

0.333048 

0.357749 

0.334287 

0.358105 

0.335228 

0.358582 

0.335861 

0.359182 

0.335932 

0.359909 

0.362126 

0.362877 

0.363635 

0.364399 

0.365170 

0.365960 

0.368500 

0.398694 

0.369249 

0.398800 

0.369862 

0.398959 

0.370329 

0.399173 

0.370644 

0.399441 

0.370677 

0.399764 

0.400746 

0.401077 

0.401410 

0.401744 

0.402081 

0.402419 

0.403493 

0.418763 

0.403815 

0.418826 

0.404079 

0.418920 

0.404282 

0.419046 

0.404419 

0.419204 

0.404433 

0.419394 

0.419971 

0.420166 

0.420362 

0.420558 

0.420756 

0.420955 

0.421576 

0.428368 

0.421762 

0.428414 

0.421916 

0.428484 

0.422034 

0.428578 

0.422114 

0.428696 

0.422122 

0.428837 

0.429268 

0.429413 

0.429559 

0.429706 

0.429853 

0.430002 

0.430462 

0.441654 

0.430600 

0.441683 

0.430714 

0.441726 

0.430802 

0.441784 

0.430861 

0.441857 

0.430867 

0.441945 

0.442211 

0.442301 

0.442392 

0.442483 

0.442575 

0.442667 

0.442951 

0.443036 

0.443106 

0.443160 

0.443197 

0.443201 
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0.455686 

0.455687 

0.455696 

0.455711 

0.455730 

0.455754 

0.455783 

0.455812 

0.455841 

0.455871 

0.455901 

0.455931 

0.455961 

0.455992 

0.456023 

0.456054 

0.456085 

0.456117 

0.456145 

0.456168 

0.456186 

0.456198 

0.456199 

0.000000 

0.003130 

0.013824 

0.029582 

0.050062 

0.074802 

0.103219 

0.130240 

0.155860 

0.179937 

0.202284 

0.222896 

0.241772 

0.258919 

0.274350 

0.288078 

0.300126 

0.310517 

0.317930 

0.323057 

0.326440 

0.328478 

0.328755 

0.003128 

0.005248 

0.014776 

0.028894 

0.047394 

0.069991 

0.096331 

0.121816 

0.146442 

0.170204 

0.192927 

0.214376 

0.234537 

0.253398 

0.270953 

0.287196 

0.302128 

0.315751 

0.326106 

0.333741 

0.339090 

0.342478 

0.342908 

0.006040 

0.119448 

0.007780 

0.016799 

0.030195 

0.047807 

0.069419 

0.094760 

0.143479 

0.166850 

0.189543 

0.211273 

0.231916 

0.251457 

0.269884 

0.287186 

0.303358 

0.318393 

0.330054 

0.338813 

0.345050 

0.349050 

0.349537 

0.007812 

0.009253 

0.017833 

0.030604 

0.047445 

0.068195 

0.092652 

0.116621 

0.140101 

0.163087 

0.185579 

0.207502 

0.228555 

0.248690 

0.267892 

0.286148 

0.303448 

0.319781 

0.332650 

0.342455 
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